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I. INTRODUCTION 

Following the discovery of the top quark at the Tevatron in Run I , one of the aims 

of the current round of data-taking is to study the top quark in more detail. In addition to 
the accumulation of further statistics in the tt pair-production channel, both collaborations 
are perform a S e arcll for single top product ffflflflH. Since S ,n gl e top proton 
proceeds by the exchange (or production) of a W boson, it offers another window into the 
weak interactions of the top quark and potentially can lead to a direct measurement of Vtb- 
Relative to the tt pair-production channel, single top production is suppressed by the weak 
coupling but favored by phase space. 

We shall report here on the method of inclusion of s ingl e top processes into the general 
next-to- leading order Monte Carlo program M CFM ^, Q, 11 1 and give phenomenological 



results relevant to the search strategy in Run II. At the Born level the single top processes 



which we include are the s-channel process 



14] 



u + d -> W* -> t + b, 
the t-channel process [l^,[l^,[lS, 12, 3], 



(1) 



b + u — > t + d, 



(2) 



and the tW mode 



19 



20] 



b + g ^t + W~. (3) 

At the Tevatron, the tW mode, Eq. Q, evaluated in the Born approximation represents 
less than one percent of the total single top cross section; in this paper we will not consider 
it further. The processes shown in Eqs. are schematic. The actual implementation in 
the program includes a sum over all contributing partons in the initial state. In addition, 
we also include the leptonic decay of the top quark, 



t -> v + e + + 6, 



(4) 



which allows a better comparison with experimental studies. 

Although some information can be extracted from Born-level calculations, the first serious 
approximation in QCD is obtained by including 0{pts) radiative corrections. We shall refer 
to this as next-to- leading order, (NLO). It is only in NLO that a calculation gives any 



information about the choice of factorization and renormalization scale. Moreover, if the 
calculation includes jets in the final state, it is only in NLO that one obtains first information 
about the structure of the jets. The NLO corrections to the inclusive s-channel mode 



have been presented in ref. 
been considered in ref. 1 22. 



2 If and the corrections to the inclusive t-channel process have 
I. The NLO corrections to the differential distributions for 



the production of single top, (without the decay of the top quark), which are needed for 



24j for the processes 



comparison with experimental results, have first been considered in ref. 
in Eqs. (0) and (0). For a recent update of this work see ref. j^ . 

In this paper we extend this program by adding the leptonic decay of the top quark with 
full spin correlations, as noted above, and also by including the effects of gluon radiation in 
the decay. The approximations employed in incorporating the QCD corrections are described 
in Section ITT1 

Section ITTT1 outlines the calculation of the 0(a s ) corrections to the decay of a free top 
quark which is helpful to establish notation. The subtraction method which we use to include 
the radiative corrections to the decay of a top quark in the single top production processes 
is described in Section HVl 

The search for single top is expected to be more challenging than top-antitop associated 
production, both because of the smaller cross section and the presence of larger backgrounds. 
We therefore give a full account of the signal and background processes in Section |V] Both 
the signal and the dominant background processes are evaluated at next-to-leading order. 



II. QCD CORRECTIONS 

In order to describe the inclusion of radiative corrections we shall discuss the s-channel 
process. The diagrams for the t-channel process can easily be constructed by crossing. They 
are obtained by reading the diagrams shown in Figs. from bottom to top, instead of 
from left to right. Some of the statements made below are specific to the s-channel process, 
but the extensions to the t-channel process are obvious. 

We shall work in the on-shell approximation for the top quark. Thus every diagram 
considered has one top quark exactly on its mass shell. Diagrams without an on-shell top 
quark are suppressed by T t / m t where T t and m t are the width and mass of the top quark. In 
this approximation the real radiative corrections fall into three types: radiation associated 
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FIG. 1: Initial-state radiation in the production stage. 



b b 




FIG. 2: Final-state radiation in the production stage. 

with the initial state (Fig. [TJ, radiation in the final state associated with the production of 
the top quark (Fig. |2J), and radiation in the final state associated with the decay of the top 
quark (Fig. EJ). The double bars in these figures indicate which top quark is on its mass 
shell. By producing the top quark strictly on its mass shell we are assured that the diagrams 
in Figs. 121 and El are separately gauge invariant. 

In addition to these real radiation diagrams we also have virtual radiation diagrams. In 
our approximation these again fall into the same three categories; virtual radiation in the 
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FIG. 3: Final-state radiation in the decay stage. 



b 




FIG. 4: Virtual radiation in the initial state. 



initial state Fig. final-state virtual radiation in the production stage Fig. El and final-state 
virtual radiation in the decay stage, Fig. El Self energy contributions on massless external 
lines have not been displayed. They give rise to scaleless integrals, which are set equal to 
zero in the dimensional regularization scheme. 

We have neglected the interference between real radiation in production and decay di- 
agrams. An example is shown in Fig. We also neglect the virtual radiation diagrams 
that link the production and decay stages, Fig. |HJ The physical reason for the neglect of 



these diagrams has been provided in Refs. j, |28( . The characteristic time scale for the 



production of the tb pair is of order l/m t while the time for the decay is 1/T t . Therefore 
in general, radiation in the production and decay stages are separated by a large time and 
the interference effects average to zero. The potentially dangerous region for this argument 
is the one in which the emitted gluon is soft since this effect is not confined to a time of 
order l/m t . The phase space for soft radiation is limited by the region in which the prop- 
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FIG. 5: Final-state virtual radiation in the production stage. 




FIG. 6: Final-state virtual radiation in the decay stage. 

agators remain resonant and the gluon energy is less than of order r t . However because of 
the cancellation of real and virtual radiation the region of soft radiation is not especially 
privileged. Thus for infrared safe variables these interference effects are expected to be of 
order a s T t / m t- 

Confirmation of this suppression for the s-channel process is provided by the work of 
Pittau who included interference between the decay and production stages. The final 
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FIG. 7: Interference between radiation in production and decay, real terms. 




FIG. 8: Interference between radiation in production and decay, virtual terms 
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results are consistent with an effect of order T/m t . A similar study has been performed 
for e + e~ — > ti including the subsequent decay of the top quarks 3(|. Here the effect of 
nonfactorizable corrections in the invariant mass distribution of the top was found to be 
very small. 

The implementation of the cancellation of soft and collinear radiation contributions which 



are separately divergent is performed using the subtraction method 



. In our program 



MCFM, we have consistently used the 



as developed by Catani and Seymour 



dipole subtraction method for massless particles 



32 



For the case of single top production we 
have a massive quark in the final state, so we have implemented a generalization of this 
scheme as suggested in j^. A useful further generalization has been suggested by Nagy 



and Trocsanyi 



35 



who introduced a tunable parameter a which controls the size of 
the subtraction region. We have extended the massive results of ref. [34] to include this 
parameter. Further details may be found in Appendix [X] 

In order to deal with radiation in the decay stage of the process we have developed a 
specialized subtraction procedure, which we discuss in Section IIVI We expect that this 
method may be useful in other contexts. In particular, it could be applied to the decay of 
the top quarks in the ti production process. 



III. RADIATIVE CORRECTIONS TO TOP DECAY 

One of the new results in this paper is the inclusion of the QCD corrections to the decay 
of a top quark in a Monte Carlo program. Since the phase space for the decay of an on- 
shell top quark factorizes from the production phase space, most of the features of the full 
calculation are present in the decay of an isolated top quark. For simplicity, and in order 
to introduce the notation, we first reproduce the well-known result for radiation from a free 
top quark j^. 

The lowest order process is t(p t ) — > W(pw) + b(Pb), where the momenta carried by the 
fields are shown in brackets. The matrix element summed and averaged over initial spin and 
color is, 

^ | 2 = 2^(l-r 2 )(l + 2r 2 ). (5) 

where we have defined r 2 = (p t — pb) 2 /mf. For the case of an on-shell W boson, r 2 = 
Myy/mf. Here and throughout this paper, the mass of the 6-quark is set equal to zero. The 
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FIG. 9: Virtual radiation in free top decay 



b 



corresponding Born-approximation width is given by. 

To 



r 2 ) 2 (l + 2r 2 ) 



(6) 



A. Virtual corrections 



The form-factor for the process t — > W + b including virtual gluon corrections is defined 

by, 

^(p b ,Pt) = u( Pb )F^ Pb , Pt )u( Pt ) (7) 

with a massless b quark and a massive t quark and momentum transfer P yy = Pt — Pb . The 
contributing diagrams are shown in Fig. |H1 The result for the form factor evaluated through 



order a„ is well known 



24, 



a s C F 



1 + 



AlTfl 2 



C 



+ 



47T 



47rr(l - e) \m 2 t - is 

+ 0(e 



n Pb 1/0 Pw 



m t 



where 



Co 



1 1/5 , . 9x 

_ ( 21n(1 _^ 



2 Li 2 r 

2 6 V y 



+ 3 m(l - r 2 ) - 2 ln 2 (l - r 2 ) - — h(l - r 2 )) 



ln(l -r 2 



(8) 



(9) 



and 

7r/l = ^(1±75). (10) 

The C*2 term will not contribute to physical amplitudes. The ultra-violet and infra-red 
divergences have been regulated by continuing to d = 4 — 2e dimensions. The dimensional 
regularization scheme is determined by the parameter r\ as follows, 

77 = 1, 't Hooft-Veltman scheme 0] 

rj = 0, four-dimensional helicity scheme . (11) 
The final result for the virtual correction to the total decay width is, 

r _. ro ^(^y( Co+ i Gi ^i) +0(£) . (12) 

We now make three parenthetic remarks which are useful for including the virtual correc- 
tions in the production processes, Eqs. (1 1121) . First, we note that Eq. © is valid for r 2 < 1. 
For the case when r 2 > 1, which is needed for the virtual corrections to the s-channel 
production process, we apply the transformation, 

7T 2 

Li 2 (r 2 ) = Li 2 (l - r 2 ) - ln(l - r 2 ) lnr 2 (13) 

6 

and analytically continue by replacing ln(l — r 2 ) with ln(l — r 2 — is). Second, we note that 



the result for the light quark vertex inc 
for the case of production, is as follows 



r M (pd,p u ) = u{p d )YiLu{p u ) 



uding the virtual corrections, which is also needed 



(14) 



a s C F ( 47r/i 2 yr 2 3 , n ,\\ 



with q = p u — pj. Third, for the case of the production processes the singularities in these 
virtual corrections are canceled by equal and opposite singularities in the integrated dipoles. 
With a little work, using results for the integrated dipole terms in Appendix El one can 
demonstrate this cancellation. 



B. Real corrections 



The Feynman diagrams for the process 

t(p t ) -> b{p b ) + W(p w ) + g{p g ) , 
10 



(15) 



Pt Pt-Pg Pb Pt Pb+Pg Pb 

FIG. 10: Real radiation in free top decay 



are shown in Fig. El The result for the matrix element squared in four dimensions is given 
as follows 



G 



F 2 



m 2 (l -r 2 )(l + 2r 2 )(2 



Pt-Pb 



Pb-Pg Pt-Pg {.Pt-Pgf 



+ 8 r 2 + 2 {Pt - P9 ~ Pb - F9J (l + 2r 2 ) 



Applying the substitutions 



Pb-Pg Pt-Pg 



Pb-Pg = -y(l - r) y 

2 

PfPg = "^-{l - r 2 ){\ - z) 



(16) 



(17) 



the four-dimensional matrix element becomes, 

1 / 2 



Y,\M\ 2 = g 2 C F \ [_L_( 



+ 8m 2 — 



Lp fe .p g U - z 
V ( 3 



1 - z 



(Pt-P ff ) s 



2 5\ l + 2r 2 

- r 2 + 2r - - J + ^ * 

V^Ll-zV l + r 2/ 



1 



(18) 



(1-z) 

In the rest frame of the top we obtain the two particle phase space for the decay t — » W+b 



d^(Pw,Pb;Pt) = (2 t")n-i (^r5 n ( Pt -p w -p b )5(pl)5(p* w -r*m* t ) 



(4vr 



i2e 



8 (2tt) 2 (m 2 )< 

The corresponding result in the rest frame of the top for the decay t^W + b + g is 

m V), \ (l-r 2 )(l-r) 2 - 4 - K)^(47r)^ 
d$ l \pw,Pb,Pg\Pt) = — — ; d 



(19) 



32 (2tt) 4 r(i - e ; 

1 /* 2/771 a x 

dz (r 2 + z(l - r 2 ))- e 
o jo 



dy y € {y ma x - y)' 



(20) 
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TABLE I: Table of integrals denned by Eq. (1231) . 



1 2 



,y (1-2) 



4 - - ln(l - r 2 ) + 2 ln 2 (l - r 2 ) + 2 Li 2 (l - r 2 ) - 5|± + 0( e ) 



y I ■' ' - "' l ' ' ' I l-r-) 



i - 3 + 2 ln(l - r 2 ) + 7^ lnr 2 + 0(e 



■^-feSi+lnd-^-^lnr 2 



-2 1 

(l+r) 2 (1-z) 2 



I + 2 - ^ lnr 2 - 21n(l - r 2 ) + O(e) 



1 2/ \ ^ 2 PV) ]nr 2 I (l+5r 2 ) n / x 

((i-r) 2 ) (r^l) +T ^i nr 2 + 0(e) 



where the upper integration limit y ma x is given by 

_ (l + r)Ml~^) (21) 
(z + r 2 {l-z)J 

Taking the ratio of Eqs. (}20j) and ()19)1 we obtain the factorized form, 
d$®(pw,Pb,P g ;Pt) = d$ (2) (p w ,p b ;p t ) x 



(1-r) 2 (m 2 ) 1 "^) 6 



16tt 2 T(l-e 



/ -Lr\2e / * / Umax 

x (r37J y dz(r 2 + z(l-r 2 ))-* J dy y~ e (y max - y)'* 

(22) 



Values of the reduced integrals defined by 

dz {r 2 + z{l - r 2 ))- e 
Jo 



1 \ /1_|_ r \2e/' i rymax 

(j(y,z)) = {^) / dz (r 2 + z(l - r 2 ))-* dyy-*(y max -y)-*f(z,y) (23) 



are given in Table U] 

Using these integrals we can calculate the contribution to the total decay width from the 
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real diagrams. The result is, 

5tt 2 



IV,., lff> a f F -^Vfl + l^-21n(l-r 2 ))- — + 2Li 2 (l-r 2 ) 
2tiT(1 - e) V ™t / U 2 e \2 V ') 6 V ; 

r ,/ 1 2n m 2/ n 2n 2r 2 (l + r 2 ) (1 - 2r 2 ) , 2 2(7r 4 -5r 2 -4) 77 

-51n(l - r 2 ) + 21n 2 (l - r 2 ) — ^ — ^ ^lnr 2 - — J- + - 

V ' K ' (1 -r 2 ) 2 (l + 2r 2 ) (1 + 2r 2 )(l - r 2 ) 2 

(24) 

The 77-dependence has been restored in this equation, so that Eq. (|2*4*j) is valid also for the 
't Hooft Veltman scheme. 



C. Result for the corrected width 

Using the virtual corrections of Eq. (jl2j) and the real corrections of Eq. (|2*4*|) we can 
calculate the value of the top width at 0(ag). We write the correction to the decay rate in 
the form 

r = r + a s T 1 (25) 
Our result is in agreement with the m 2 /m 2 limit of the original calculation in ref. 



a s ri a s 

— — — — Op 

r n 2tt 



-tt 2 + 4 Li 2 (r 2 )---2 ln(— — -) +2 lnr 2 ln(l-r 2 
3 2 V 1 — r 2 / 

(22 -34r 2 ), , (3 + 271n(l - r 2 ) -41nr 2 ) 
lnr 2 + 



(26) 



3(1 -r 2 ) 9(1 -r 2 ) 2 9(1 + 2r 2 ) 

where r = My^jm t . For r ~ 4/9 the QCD correction amounts to 

^ « -0.8a s (27) 
t 

which lowers the leading order result for the width by about 10%. 

IV. FACTORIZATION OF SINGULARITIES IN TOP QUARK DECAY 

We wish to construct a counter-term for the process 

t -> W + b + g (28) 

which has the same soft and collinear singularities as the full matrix element. This counter- 
term takes the form of a lowest order matrix element multiplied by a function D which 
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describes the emission of soft or collinear radiation, 

\M(. . .p t ,Pw,Pb,Pg)\ 2 -> \M (. - .Pt,Pw,Pb)\ 2 x D(p t .p g ,p b .p g} m 2 ,M w ) , (29) 

In the region of soft emission, or in the region where the momenta p g and pb are collinear, 
the right hand side of Eq. (}29j) has the same singularity structure as the full matrix element. 
The lowest order matrix element Aio in Eq. (|29jl is evaluated for values of the momenta 
Pw and pb modified to absorb the four-momentum carried away by the gluon, and subject 
to the momentum conservation constraint, pt — > pw + Pb- The modified momenta denoted 
by a tilde are also subject to the mass-shell constraints, pi = and p\y = p^. The latter 
condition is necessary in order that the rapidly varying Breit-Wigner function for the W is 
evaluated at the same kinematic point in the counterterm and in the full matrix element. 
We define pw by a Lorentz transformation, p^y = A^pyy fixed in terms of the momenta pw 
and p t . Because pw and pw are related by a Lorentz transformation the phase space for the 
subsequent decay of the W is unchanged. 

The general form of a Lorentz transformation in the plane of the vectors pt and pw is 
given by (p£, = A^p w ), 

a uv sinh(x) / \ 

A ^ = 9»v + . ===== [ptPw ~ PwPt ) 

V(Pt -PwY -p w Pt 

+ t C ° Sh ^ — (ft ■ p w (PtPw + PwPt) - Pw PtPt - Pt PwPw) ( 30 ) 
{Pt ' Pw) — PwPt v ' 

The transformed momentum of the b quark is fixed by pb = Pt~ Pw- F° r the special case in 

which we impose the condition pi = (p t — pw) 2 = we get 

1 



sinh(x) = g 2 - (p^ - p 2 w )pt ■ Pw + {p 2 + Pw) J (Pt • Pw) 2 ~ PwPt 
A PtPw L 

cosh(x) = I 2 + (p 2 + p 2 w )pt ■ Pw - {p 2 - P 2 w)\l (Pt ■ Pw) 2 ~ PwPt ( 31 ) 
Acting on the vector pw the Lorentz transformation becomes 

PW = OL (p w - — ^ Pt) + P Pt ( 32 ) 
Pt 



where the constants are given by 



a = - P * P " w = (33) 
2 V (Pt ■ Pw) 2 ~ PwPt 

P = P l±j^ (34) 
14 



Eq. (|32|) makes it clear that, in the top rest frame, the transformation on pw is a Lorentz 
boost along the direction of the W. 

The momenta of the decay products of the W can similarly be obtained by the same 
Lorentz transformation, Eqs. (|3UI31j) . 



A. Subtraction counterterm 



From Eq. (|2*2*|) we may write the phase space for the decay of an on-shell top quark as 



d® (3) (pw,Pb,P g ;Pt) = d^ 2) (p w ,p b ;p t ) I [dg(pt,pw,y,z)) 

= d$ {2) (p w ,p b ;p t ) J[dg(p t ,p w ,y,z)} (35) 

The equivalence in Eq. ()35jl follows from Eq. because since pw and 

Pw are related by a boost. From Eq. (j22j) we see that the phase space integral for the emitted 
gluon is given by, 



[dg(PuPw,y,z)] 



(l-O 2 , 2 Ul-e) (47T) 



16vr 2 v *' 



1 — r 



r(i-e) 

,2 



/ ^(r 2 + z(l-r 2 ))- / ^y- e (z/ m a,-y)- £ (36) 
Jo Jo 



where y and z are given by Eq. ()17|) . By extension of Eq. (j!8|) we choose the counterterm 
to be, 



D((p t + Pg y, (pi+PgY, mf, M 2 W ) = g 2 ^C F 



Pb-Pg\{l ~ Z) 



-l — z—r]e(l — z] 



it.gf 



(37) 



where the role of the parameter 77 is defined in Eq. fill}) . Performing the integral using the 
results of Table |1J we obtain the following result for the integrated counterterm, 



[dg(pt,Pw,y,z)} x D((p t +p g ) 2 ,(p b + pg) 2 ,m 2 t ,Mw] 



asCp (4:71/1 



2\c 



2tt mf r(l - e) 
1 



i + -f|-21n(l-r 2 



e V2 



25 1 

+ T + 2~V(l-r 2 ) 2 (1-r 2 
5tt 2 



+ 7 In r 2 



2 1 -r 2 



+ 2Li 2 (l -r 2 ) 



?7 



6 



51n(l-r 2 )+21n 2 (l-r 2 ) + - (38) 
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TABLE II: LO and NLO cross sections for single top-quark production at the Tevatron and LHC for 
mt = 178 GeV. The branching ratio for the decay of the top quark is not included. Cross sections 
are evaluated with CTEQ6L1 (a s (M z ) = 0.130) and CTEQ6M (a s (M z ) = 0.118) PDFs 48], and 
all scales set to mt- The errors represent Monte Carlo statistics only. 



Process 


[TeV] 


olo [pb] 


(Jnlo [pb] 


s-channel, pp (t) 


1.96 


0.270 


0.405 ±0.0003 


s-channel, pp (t) 


14 


4.26 


6.06 ±0.004 


s-channel, pp (t) 


14 


2.58 


3.76 ±0.003 


i-channel, pp (t) 


1.96 


0.826 


0.924 ±0.001 


i-channel, pp (t) 


14 


146.2 


150.0 ±0.2 


i-channel, pp (t) 


14 


84.8 


88.5 ±0.1 



TABLE III: LO and NLO cross sections for single top-quark production at the Tevatron and LHC 
for mt = 178 GeV. The branching ratio for the decay of the top quark is not included. Cross 
sections are evaluated with the MRST2002 NLO PDF set with a s (M z ) = 0.1197, and all 
scales set to mt- The errors represent Monte Carlo statistics only. 



Process 


[TeV] 


olo [pb] 


(JNLO [pb] 


s-channel, pp (t) 


1.96 


0.285 


0.404 ±0.0003 


s-channel, pp (t) 


14 


4.57 


6.17 ±0.004 


s-channel, pp (t) 


14 


2.85 


3.86 ±0.003 


i-channel, pp (t) 


1.96 


1.009 


1.032 ±0.001 


i-channel, pp (t) 


14 


160.1 


154.4 ±0.2 


i-channel, pp (t) 


14 


96.9 


92.4 ±0.1 



V. PHENOMENOLOGICAL STUDIES 



The method of the previous section has been implemented in the Monte Carlo program, 
MCFM, allowing us to make predictions for kinematic distributions of both signal and 
background events. Before proceeding to describe our results for jets and the decay products 
of the top, we present a number of results on total cross sections. 



A. Total cross section 



Tables ILTl and IIHI give total cross sections using recent parton distributions and the up- 
dated top quark mass 4|| . The input parameters which we use throughout this phenomeno- 
logical section are presented in Table IIVI We have not performed a full analysis of the 
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TABLE IV: Input parameters 



m t = 178. GeV T = 1.651 GeV 

M w = 80.4 GeV T w = 2.06 GeV 

M z = 91.188 GeV T z = 2.49 GeV 

G F = 1.16639 x 10~ 5 GeV" 2 a' 1 = 132.351 

g^y = 0.42651 e 2 = 0.0949475 



theoretical errors on these predictions. It is important to remember that the results for 
the t-channel process depend on the 6-quark parton distribution, which is calculated rather 
than measured. This is the source of the relatively large difference between the t-channel 
Tevatron cross sections in Tables |H] and 11111 compared to the s-channel process which is 
dominated by well-constrained valence quark distributions. 

We have also compared our leading order (LO) and next-to-leading order (NLO) results 
for the total cross sections (without radiation in the decay) at y/s = 1.96 and 14 TeV with 



the results in ref. 24J. With the appropriate modification of the input parameters, we find 
agreement with their results within the quoted errors. 

We now investigate what effect the inclusion of QCD corrections in the decay of the top 
quark has on the total cross section. Radiation in the decay should not change the total 
cross section. However in a perturbative approach this means that the difference should be 
of higher order in a s . Including QCD radiation effects only in the production, we obtain 

aB t ^ hve = aB w ^ ve = a^)B w ^ ue + a s a(i)B w ^ ue (39) 

where ovo) and <T(i) are the LO and NLO cross sections, because in this order the width 
t — > bW is equal to the total width. When we also include radiative corrections to the 
decay, we use the radiatively corrected total width T = T + a s Ti, cf. Eqs. (I6|25jl . 



a Bt^bue+x = cr ( ) h a s a (i) h a ( ) (4U ) 

The difference between these two expressions is of 0(a 2 ) and is given by 
D d T B w ^ ue 

V-Dt-tbve+X — V&t-^bve = CK B (T(i) a s (J(i)nw^ve 

= -a^W f ( 41 ) 

As shown in Table IVl the numerical differences are less than 3% for the s-channel process 
and less than 0.5% for the t-channel process. 
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TABLE V: Comparison of LO and NLO cross sections for single i-quark production at the Tevatron 
and LHC. The NLO calculation is performed both without including QCD effects in the decay 
(aBt-^bue) and also when it is included (aB t ^bue+x)- The top quark mass is mt = 178 GeV and 
cross sections are evaluated using MRST2002 NLO PDFs with all scales set to mj. The errors 
represent Monte Carlo statistics only. Note that the values of Tt at LO and NLO are 1.6511 GeV 
and 1.5077 GeV respectively and the branching ratio of the W into leptons is Br(VF — * ev) = 0.1104. 
The differences in the rates are in good agreement with the formulae in Eq. (|41j) . 



Process 


[TeV] 


a B t ^ bue (fb) 


crBt^bue (fb) 


aBt^bve+x (fb) 


s-channel, pp (t) 


1.96 


31.54 ±0.01 


44.64 ±0.03 


45.88 ±0.03 


s-channel, pp (t) 


14 


503.9 ±0.2 


681.0 ±0.4 


698.7 ±0.4 


i-channel, pp (t) 


1.96 


111.34 ±0.06 


113.95 ±0.12 


113.96 ±0.12 


i-channel, pp (t) 


14 


17690. ±8 


17048. ±16 


16975. ±16. 



B. Signals and Backgrounds at NLO 

We shall consider the signal for single top production to be the presence of a lepton, 
missing energy and two jets, one of which is tagged as a 6-jet. In the case where we have 
two tagged jets, we choose the jet to be assigned to the top quark at random. For clarity, 
we shall describe the processes at the parton level choosing specific partons in the initial 
state. In our program we sum over all the species of partons present in the initial proton 
and antiproton. For reactions considered at NLO, there can also be additional partons in 
the final state. We use the Run II /^-clustering algorithm to find jets, with a pseudo-cone 
of size R = 1.0 (47]]. The first reactions to consider are the two signal processes 

u + d — > t + b 

I (42) 
U v + e + + b 

and 

b + u — > t + d 

I • (43) 

U v + e + + b 

Note that we present numerical results for the processes shown in Eqs. (|42I43|) . summed 
over species of initial partons, i.e. the production of a t-quark (rather than a t) with the 
decay t — > u + e + + b + X. Thus in a hypothetical experiment with equal perfect acceptances 
for electrons of both charges and for muons of both charges the signal (and background) 
will be four times bigger, (if we ignore possible signatures coming from W — > rv T \ Unless 
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explicitly stated otherwise, our NLO results for the signal processes are calculated including 
QCD corrections in both the production and decay of the top quark. 

The background processes which we consider are of several types. The irreducible back- 
grounds are, 

u + d -> W + + b + b 
U v + e + 



(44) 



u + b -> W + + d + b 
U v + e+ 



u + d -> W + + Z 



z/ + e" 



6 + 6 



(45) 



(46) 



There are also backgrounds related to tt production which contribute to the W+ 2 jets 
process. The case where both top quarks decay leptonically, 



(47) 



u + u — > t + t 

U b + e~ + v 
v + e + + b 

contributes if the electron, e~ (or muon, fj,~) fails the cuts. If one of the top quarks instead 
decays hadronically then there can also be a contribution when only two jets are observed, 
either because of merging or because the extra jets lie outside the acceptance (or both), 

u + u — > t + t 

U b + q + q (48) 
> v + e + + b 

A significant background process involves W + two light jet production where one of the 
light quark jets fakes a 6-quark, 

u + d W + + 2 jets 
U v + e+ 

Further backgrounds involve the mistagging of a c-quark as a 6-quark, 

u + s — > W + + u + c 
U v + e + 



(49) 



(50) 
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TABLE VI: Cuts for the single top analysis presented in this section. 



Lepton pt 


p e T > 20 GeV 


Lepton pseudorapidity 


\r] e \ < 1.1 


Missing Ex 


E T > 20 GeV 


Jet pt 


^ T et > 15 GeV 


Jet pseudorapidity 


|^ et | < 2.8 


Mass oib + l + v 


140 < m b iu < 210 GeV 



u + d -> W + + c + c 

I (51) 
U i/ + e + 

I (52) 
U z/ + e + 

Our estimation of the rates for these processes depends on the cuts shown in Table IVI| 
which have been chosen to mimic those used in an actual Run II analysis |4j. As well as 
the normal jet and lepton cuts, we perform a cut on the missing transverse energy and the 
mass of the '& + I + ^'-system, muv The missing transverse energy vector is the negative 
of the vector sum of the transverse energy of the observed jets and leptons. The mass of 
the putative top system is determined by reconstructing the W and combining it with the 
tagged 6-jet. In reconstructing the W the longitudinal momentum of the neutrino is fixed 
by constraining the mass of the ev system to be equal to Myy. The two solutions for the 
longitudinal momentum of the neutrino are, 
1 



Pi 



pi {Ml - Mf + 2\p T \\p T \) ± E e J(M$y - M T )(Myy - m\ + ±\p T \\p T \) (53) 



2|^| 2 

In this equation p T is the measured missing transverse energy and 

M 2 T = (\p T \ + \p T \f - (pS + Vt? (54) 

is the transverse mass of the W. We resolve the twofold ambiguity in p u L by choosing the 
solution which gives the largest (smallest) neutrino rapidity for the W + {W~). 

Our results are shown in Table rvTH In this table we first show the cross sections calculated 
assuming a perfectly efficient detector, both at LO and NLO where possible. The column 
labeled 'Efficiency' gives the rescaling factor that should be applied in order to take account 
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TABLE VII: Cross sections for the Tevatron Run II {jpp, \fs = 1.96 TeV) in femtobarns, calculated 
with MRST2002 and using the renormalization and factorization scale [i. The cross sections contain 
two jets, subject to the cuts of Table IVT1 Next-to-leading order cross sections are also shown, where 
known. The last column estimates the importance of signal and background processes by including 
nominal tagging and mistagging efficiencies, using the NLO cross section where possible. 



Process 


Scale 


CTLO [fb] 


0JVLO [fb] 


Efficiency 


<7[fb] 


s-channel single top, Eq. (fl2*|) 




10.3 


11.7 


i _ 

J. 


n - f^ 2 


7.4 


s-channel (with decay radiation) 


m t 


10.3 


11.3 


1 - 


(l-^) 2 


7.2 


i-channel single top, Eq. (|4*3*)) 


mi 


38.8 


29.4 


£b 




11.8 


i-channel (with decay radiation) 


m t 


38.8 


26.6 


£b 




10.6 


Wbb, Eq. (05) 


M W 


36.0 


47.5 


1 - 


(1 - e b ? 


30.4 


W + bj, Eq. @HJ) 


M w /4 


26.5 




£b 




10.6 


WZ, Eq. (H)J) 


M w 


3.64 


3.91 


1 - 


(l~e b ) 2 


2.5 


ti, Eq. (gZJ) 


m t 


4.34 




2 x 


(1-(1-^) 2 ) 


5.6 


tt, Eq. (@E1) 


m t 


4.94 




1 - 




3.2 


W+2 jet, Eq. (gHJ) 


M w 


5530. 


7030. 


1 - 


(1 " h? 


35.1 


us -> Wuc, Eq. (EOJ) 


M w 


324. 




fc 




19.4 


Wcc, Eq. (EU 


M w 


36.0 


47.5 


1 - 


(1 " fc? 


5.5 


W + cj, Eq. © 


M w /A 


54.7 




/e 




3.3 



TABLE VIII: Nominal efficiencies used throughout this section. 


6-tagging efficiency 


£b = 


40% 


c-mistagged as b 


fc = 


6% 


Jet fake rate 


fj = 


0.25% 



of experimental efficiencies and fake rates. To give some idea of the importance of the 
background processes, we have used the nominal (and perhaps optimistic) values of these 
quantities shown in Table IVIIII to obtain the final column of Table IVHI 

A strategy for searching for the single top processes at the Tevatron involves examining 
the H T distribution 4], where H T is denned by, 

H T = |p T (lepton)| + \p T \ + |Pr(jet)| . (55) 

In Fig. ^2 we show the distributions of the signal and background processes described 
above, rescaled by the nominal efficiencies given in Table jVITTI As in the table, the NLO 
calculation is used wherever it is known and the signal processes include gluon radiation 
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FIG. 11: Predictions for Ht showing single top signal, (top left), irreducible backgrounds (top 
right), tt backgrounds, (bottom left) and jet and charm fake rate (bottom right). Next-to-leading 
order estimates are used for the cross sections where available. 



in the decay of the top quark. We see that the Ht distribution is harder for the single 
top distributions, than all backgrounds, except those which involve top production (which 
are however, small). Despite this fact, one sees that the background rates are still large 
in the region of Ht where the signal processes peak. This is demonstrated more clearly 
in Figure ^1 where we show the Ht distribution for the sum of all the backgrounds, the 
sum of the two single top processes and the total when combining signal and background. 
Although the single top signal represents about 50% of the events in the bins of interest, an 
observation of this top production mechanism using the Ht distribution is heavily reliant 
on accurate predictions of both the rates and shapes of the background processes. We stress 
that Figs. ^2 and ^] depend on the particular values chosen for the efficiencies and will 
improve if better values are achieved. 
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H T [GeV] 

FIG. 12: The Ht distributions of signal, background and signal plus background. The curves 
correspond to the sum of the appropriate distributions shown in Figure 1111 



A further quantity that may be used to discriminate between the t-channel signal process 
and backgrounds is Qrj, where Q is the charge of the lepton (in units of the positron charge) 
and 7] is the pseudorapidity of the untagged jet. In this analysis |j| an additional cut is 
applied, requiring that one of the two jets must have a px > 30 GeV. This extra cut serves 
to reduce the dominant backgrounds in Table IV 111 by 15-30% whilst leaving the signal 
virtually unchanged. Since this is a targeted search for the t-channel process we have also 
rejected events with two tagged jets, since at leading order, Eq. (j43|) contains only one b- 
jet. Results for the Qr] distribution are shown in Fig. demonstrating that the t-channel 
process is enhanced at large pseudorapidity compared to the backgrounds. As before, even 
the pronounced effect observed in this distribution can be hidden when comparing the signal 
to the sum of all backgrounds, as shown in Fig. 

Finally, we end this section with a comment on the effect of including the radiation in the 
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FIG. 13: Predictions for Qrj showing single top signal, (top left), irreducible backgrounds (top 
right), tt backgrounds, (bottom left) and jet and charm fake rate (bottom right). Next-to-leading 
order estimates are used for the cross sections where available. 



top quark decay at NLO. We find that although the overall rate is lowered (cf. Table l*VTT|) . 
the shapes of the H T and Qi] distributions are not altered significantly. This is demonstrated 
in Fig. EE where we compare the Ht (s-channel) and Qrj (t-channel) distributions with and 
without radiation in the decay. However this is a feature of our specific choice of cuts and 
we do not know this to be true in general. 



VI. CONCLUSIONS 



We have presented the first results of our NLO Monte Carlo program which describes 
the signal and background for single top production. The new feature of our analysis is the 
inclusion of the decay of the top quark. Since our analysis is performed at next-to-leading 
order, we have included radiative effects both in the production and in the decay. Radiation 
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FIG. 14: The Qrj distributions of signal, background and signal plus background, 
correspond to the sum of the appropriate distributions shown in Figure 



The curves 



in the decay is performed using a new subtraction method, which may be useful in other 
contexts. 

We have studied the effects of the NLO corrections on two distributions that are relevant 
for single top searches at the Tevatron. The inclusion of the top quark decay is imperative 
for such an analysis for two reasons. Firstly, the cuts that are applied (in order to match 
the experimental analysis) are on the decay products and not on the top quark itself. Sec- 
ond, one of the distributions of interest, H?, can only be calculated when the momenta 
of the decay products are known. We have also performed the calculation of a number 
of important backgrounds at NLO and examined the feasibility of using variables such as 
Ht and Qi] to discriminate between the single top signal and the main backgrounds. We 
find that the inclusion of radiation in the decay makes very little difference in the shapes 
of these distributions, but further decreases the exclusive two-jet signal cross section. Our 
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FIG. 15: Comparison of the NLO distributions for Ht in the s-channel process (left) and Qr] in 
the i-channel process (right), with and without radiation in decay. 



treatment is obviously approximate, because we assume pr-independent efficiencies, stable 
b and c quarks and include no showering or hadronization. Nevertheless it confirms that the 
search for single top is extremely challenging and suggests that significantly larger tagging 
efficiencies and/or methods with greater discriminatory power will be needed to observe sin- 
gle top production. We hope to have contributed to this search by providing more reliable 
information on the kinematic structure of both the signal and the background. 
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APPENDIX A: INTEGRATION OF DIPOLES 

As described in the text of the paper the production cross section was calculated using 
a subtraction method The particular subtraction terms were taken from ref. j^] (CS) 
for the subtraction terms with massless particles and from ref. |Q (CDST) for the dipoles 
involving massive particles in the final state. 

In this appendix we describe two minor modifications of the methods in those references. 
The first is that we allow ourselves the freedom to work in the four-dimensional helicity 
scheme. This scheme has been used for the calculation of many virtual matrix elements, so 
we believe that our results may be useful in other contexts. 

The second modification is that, following Nagy and Trocsanyi we introduce a 

tuneable parameter a which can be used to reduce the range of integration over the singular 
variable, a = 1 corresponds to integration over the whole dipole phase space and for a < 1 
the volume of phase space for the dipole subtraction is reduced. As detailed in this 
serves a number of purposes. Since the subtraction is smaller, the mis-match between real 
matrix elements and subtractions is also reduced, leading to fewer events where the event 
and counter-event fall into different bins. In addition, for small a the counter-terms are 
only calculated in events where they actually play a role in canceling a singularity. This 
can lead to a saving in computer time. Lastly, the lack of a-dependence in the results is a 
valuable check of the numerical implementation. For the massless case the integrations of 
the dipoles in reduced phase space volumes have been provided in 3(|. The results for the 
a dependence of the massive dipoles are new. 

To keep the appendix short we try and use a similar notation to refs. and refer 

back to specific equations in those references (CS Q and CDST Q) 

when appropriate. 

1. Initial-state emitter with initial state spectator 

As explained above we have used a slight generalization of the dipole phase space (CS, 
Eq. (5.151)) where the variable Vi is modified by the factor Q(a — V{). The variable Vi is the 
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rescaled value of the propagator, defined as 



PaPi 
PaPb 



(Al) 



where p a is the initial state emitter, pi is the emitted parton and p& is the other initial 
state parton which is the spectator. For a full description see section 5.5 of CS. The dipole 
integrand which we subtract is obtained by introducing r\ into CS, Eq. (5.145) for the q,q 
case and CS, Eq. (5.146) for the g, q case: 

2 



(s\V qa9i ' b {x it0b )\s') = 8ir(i 2e a s C F 5 s 



1 — x. 



i,ab 



- (1 + x i>ab ) - rje(l - x it ab)) 



(A2) 



s\V 9 ^' b (x ijab )\s') = 8vr/i 2t « s T R [1 - tje - 2x i>ab (l - x i)ab )} S s 



(A3) 



As in Eq. (jlljl for rj = 1 we are in the 't Hooft-Veltman scheme. Setting rj = we are in the 
4-dimensional helicity scheme. 

The result for the q, q case is a generalization of CS, Eq. (5.155) and is given by, 

V™(x;e,a) = C F { \ i-^ )8(l-x) 



+ 77(1 -x) - (l+x) (2 m(l - x) - - 



+ t)[l — x — a)— ^m 



[1-xl 



+ 4 



1 — x) V 1 — x 

ln(l — x 
1 — x 



0(6) 



(A4) 



The result for the g, q case (after averaging in d dimensions in the 't Hooft Veltman scheme) 
is given by, 



x 



V^(x;e,a)=T R ^(l-x) 2 + x 2 



2 ln(l - x) - - + 0(1 - x - a) In ( 

e \ 1 — x 



+ 2r)x(l-x) } +0(e) 



(A5) 



In the limit a = 1 (and 77 = 1) these functions reduce to those derived by the explicit 
expansion of CS, Eq. (5.155). They are also in agreement with the results of Nagy J43 1. 
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In the above equations we have followed the conventions of CS. However in CS the overall 
factor multiplying the functions V contains a term of the form, (cf. CS Eq. (5.152)) 



In our program we choose to write this factor as 

4:71 fi 2 \ ( 4ti[i 2 x\ ( Att/jL 2 



1 + e In x + 



(A6) 



(A7) 



ZPaPb J \2paiPb J \2PaiPb^ 

since it is the transformed momenta which are held fixed when we perform the x integration. 
In our program the additional In x terms are included in our results for the integrated dipoles. 
These additional terms containing lnx are accounted for in the paper of CS in a different 
fashion. For a similar reason, our program contains additional lnx terms in the numerical 
implementation of Eqs. (|A9|A11|) . 



2. Initial-state emitter with final-state spectator 



We have two cases which we have to consider 

a) Initial q — > q + g, with a massive final state spectator; 

b) Initial q — > q + g, with a massless final state spectator 

For ), the phase space is the generalization of CDST, Eq. (5.79) with an extra 

factor of Q(a — Zi) which implements the reduction in the phase space volume. The variable 
z + is defined by, 

I — x , , . 

z+ = -. (A8) 



The dipole integrand is given by a generalization of CDST, Eq. (5.81). The result is 

r-2 1 



I$(x;e,a) = C F { - 



L[l-x] 



— 1 — X 



1 



+ 5(l-x) 



1 7T 



6 



m(l + 



r.2\ 



+ 2 Li 2 (-/i 2 ) + 2 ln(/i 2 ) ln(l + ji 2 ) - - ln 2 (l + 



+ 4 



ln(l — x) 



y l-x 

— Q(z + — a) 

- (1 +x)\n 



,~,2\ 



ln(l + fi 

[1-^]+ "(l-a;) 



2 , /2-x + xu 2 
In 



l + jl 



r,2 



. .logl ^^ )-(l + x)ln(£±) 
l — xj \a{± — x + z+) J \ a / 



(1-a:) 2 
(1 — x + xjl 2 ] 



+ r}{l-x) }+0(e) 



(A9) 
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For a = 1 and 77 = 1 this agrees with CDST, Eqs. (5.90) and (5.88). Here we have introduced 
the variable 

A 2 = — = ipf-r- (A10) 

which only depends on p a i and pj which are defined in CDST Eq. (5.73). This is helpful 
since it is these transformed momenta which are held fixed when we do the x integration. 

For case (b), the dipole phase space is given by CS, Eq. (5.72) with an additional factor of 
Q(a — Ui). The integrand is given by the usual ^-dependent generalization of CS, Eq. (5.77). 
Performing the integration yields, 

V q ' q (x; e, a) = C F L)(1 - x) - In + + (\~ ln («) ~ ln ( X - x )) i 1 + x ) 



1-x) 



This result for rj = 1 is in agreement with CS, Eq. (5.83). It is also in agreement with the 
a ^ 1 results of Nagy [43]. Since the \i — > limit is smooth it can also be obtained from 
the limit of Eq. (jXHjl . 

3. Final-state emitter with initial-state spectator 

Once again there are two cases which we have to consider, 

a) Final Q — > Q + g, with a massless initial-state spectator; 

b) Final q —>■ q + g, with a massless initial-state spectator. 

We shall deal with the two cases in turn. The phase space for the first case is given 
by CDST, Eq. (5.48) with the addition of the factor Q(x — 1 + a) to reduce the phase- 
space volume. The integrand is given by CDST, Eq. (5.50) and yields a result which we 
choose to decompose into three separate contributions. Our decomposition of is into a 
delta-function, plus-distribution and regular parts. 

I a gQ {x\ e, a) = C F {5(1 - ar) J a J (ji Q , e, a) + J; q + (x, (jl q , a) + J^{x, fi Q , a)} + 0(e), 

where we have introduced the variable (c.f. Eq. (5.45) of CDST) 

^ = (A12) 
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This is a different choice than the one given in Eq. (5.56) of ref. and more suitable for 
implementation into our program. The three contributions to Eq. (jA12|) are 



J" Q 5 (//, e, a) 



i-iln 

e e 



r,2 
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2Li 2 (-/i 2 )- — + 2 + 21n(a) 



ln[i±^]-l 



1-2-2 , 1 i..2, 



+ -lnV 2 +2^(l + /i 2 ) - 2 ln/i 2 ln(l + /i 2 ) + ln(/2 2 ) 



A 2 



(A13) 



J a + (a;,/i,a) 



[J^(x,/i,a)] 



1 — x 



(2 - x + /i 2 ) /i 2 



+ 



In 



2(l-x + /i 2 ) 2 (1-x) \ (1 +/i 2 )(l -x + /i 2 ) 



9(x - 1 + a) 



(A14) 

where we have defined ft 2 = fi 2 /x because the x integration is performed at fixed fx 2 . The 
distribution is interpreted as follows 



dx f( x ) , , _ 
o \ l — x 



l-a 



1 ,,(/(*) -/(D) 

1-Q 1 X 



(A15) 



so that as a —>■ 1 we recover the normal plus distribution. 

Eq. (|A14j) with a = 1 should be compared with Eqs. (5.58, 5.59 and 5.60) of CDST. The 
comparison is mostly easily performed by showing that the two forms are identical for x ^ 1 
and that the integral 



dx I^ n (x:e) 



(A16) 



is identical for the two forms. 

Now we consider case b), in which both the emitter and spectator quarks are massless. 
The phase space is given by CS, Eq. (5.48) and the subtraction is given by the ^-dependent 
generalization of CS, Eq. (5.49). The result is 



+ 



V qg (x;e,a) = C F { 5(1 - x) 
2 ln(2 - x) 



1 — x 



e(x-l + a)-2 



In a) - ln 2 (a) H h 

2 2 2 2 v ; v ; 2e e 2 

ln(l - x) \ 3 / I 



1 — x I 2 \ 1 — x 

1 — a \ / 1— q. 



0(e) (A17) 



Eq. ()A17|1 is in agreement with CS, Eq. (5.57) when rj = 1. It is also in agreement with the 
a-dependent results of Nagy 
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4. Final-state emitter with final-state spectator 



We now consider final state radiation of a gluon from a quark line with a final state 
spectator. This can either be radiation off a massive line, with a massless spectator, or 
radiation off a massless line with a massive spectator: 

a) Final Q — > Q + g, with massless final-state spectator; 

b) Final q — > q + g, with massive final-state spectator. 

For case (a) we have the phase space from CDST, Eq. (5.11) multiplied by 6 (a — yij,k) 
where the and Zi integrations range over, 



< y ijt k < 1 



< Zi < 



(! - V 2 Q )yij,k 



\Pq + (! - t*Q)Vij,k] ' 
The dipole subtraction is a generalization of CDST, Eq. (5.16). 



(A18) 



(s\V gQ , q \s>) =8^a s C F \ 2 - d f- k 

I 1 - zAl - y ijt k ) Viji 



m. 



S ss > (A19) 



«jO--yij,k) v ij:k i 

Vij,k, v ij:k and defined in CDST, Eqs. (5.8), (5.14) and (5.12). Following 

CDST, Eq. (5.23) we divide the dipole subtraction into an eikonal piece and a collinear 
piece. 

I gQ , k (H, e) = C F [2/ eik (/i„ m e) + 7 coll (^, m e)] (A20) 

For a massive quark there is no collinear divergence and hence only the terms with soft 
singularities are necessary. The terms which subtract a putative mass singularity ensure 
continuity in the small mass limit. For a large enough mass they are irrelevant and they can 
be removed by setting = 0. 

The eikonal integral is defined by 

1 /4vr/i 2 ^ e 



as 

2tiT(1 - el 



Q 2 



/ eik GuQ,0;e) 



[d-Pi(Pij,Pk)] 



d>-KH 2e as 



2p iPj 1 - Zj(l - y ijtk ) 



0(a-y ijtk ), 



(A21) 
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and performing the integration we obtain, 



I^q, 0; e) = _ 2 Lia(1 _ ^ _ H ^ q) ln(1 _ ^ _ 1 ln 2 (/ ^y 



- ln(a) ln(/4) - Li 2 + Li 2 ( ^ Q 2 J ) + 0(e) . (A22) 

V \Iq / \ fig / 

Taking a = 1 we can compare with the appropriate limit of CDST, Eq. (A.l). 

The remaining terms in Eq. (|A19|) are referred to as collinear subtraction terms (even 
though the term proportional to rriq in Eq. (jA19|) has only a soft singularity). The result 
after integration is, 



r col W0;e) = i + + 2 + ln(^)(l + ^)+ln( / iyi^^- 21n (l-^ 
e z i \Iq 



+ 1 



(1-^q) w (a + (l-"K 



ln(ct + (1 — a)u%) 
- 2 ln(a) + 2 V / ^ + 0(e) . (A23) 

Setting a = 1 and <fi = 1 we recover 

J col W 0; e) = 1 + 3 + I ln(/z* ) - - 2 ln(l - ^) (A24) 

This last result can be extracted from the appropriate limit of CDST, Eq. (5.35). 

For case (b) the phase space is again taken from CDST, Eq. (5.11) but with the constraint 
Q(ay + — Hij^k) and integration limits 

1 — Mo 

< Vijjc <y+= : ■ V 
1 + Mq 

(l-v ijtk )/2< z <(l+v ij>k )/2 (A25) 



where 

y/[2^ + (l-^)(l-^ fc P -4^ 

(! -P%)( 1 -Vi3,k) 

For this case the subtraction term is 



v ,, k = vi^t^ 2\/\ yv 'T — • (A26) 



(s\V gq , Q \s') = 8n^a s C F ( f_ - ^ [1 + % + e(l - ~z 3 )r,}\ 5 S 

I %'v Vij,k) ^ij,k J 



(A27) 



The correction pieces for the a-dependent restricted range of the integral over y can be 
performed with the help of the transformation 



x = y+ - y + \ (y+ - y) (— - y) (A28) 
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as suggested in Appendix C.l of ref. 



The result for the eikonal integral is, 



/ eik (0,/i Q ;e) = — 



1 ln(l-/4) 



+ U 2 (l-fJL 2 Q ) 

1 



2e 2 



-7T 



+ ln 2 (l-^) 



12 

1 - j/i + 2x+y+ 



H ln z . 

2 V(l + ?/+-£+) (l-y+ + 2 H 



-In 2 
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+ 2 
+ Li 2 
+ Li 2 
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1 + 2/ 
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Li 2 



1 - y 2 + 2a;+i/ + 



ln(i±^lln^-^ +2 ^ + 



2^ 



Lb 



1 + J/+) 2 
1-2/ 



0(6) 



where 



-a) + J(l- a)(l - a?/ 2 



Setting a = 1 yields agreement with the appropriate limit of CDST, Eq. (A.l). 
The contribution from the collinear piece which has no soft singularity is 



/ coll (0,/i Q ;6) 



2i" 31a ( 1 -^) + — ■ i-^ 



- 2 



/^q(1 ~ g/fg) 
(1 - ^) 



3 

2 L 



ln(a)+2/+(l-a) + 0(e) 



(A29) 



(A30) 



(A31) 



After setting a = 1 and r] = 1, this can be compared to CDST, Eq. (5.35). The total 
contribution for the integrated dipole is defined in Eq. (jA20|) . 



[1] F. Abe el al. [CDF Collaboration], Phys. Rev. D 50, 2966 (1994). 

[2] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 74, 2626 (1995) arXiv:hep-ex7 9503002 . 
[3] S. Abachi et al. [DO Collaboration], Phys. Rev. Lett. 74, 2632 (1995) |arXiv:hep-ex/ 9503003 . 
[4] A. Juste [CDF Collaboration] , |arXiv:hep-ex/040604l"| 

[5] B. Abbott et al. [DO Collaboration], Phys. Rev. D 63, 031101 (2001) arXiv:hep-ex/0008024 . 

[6] V. M. Abazov et al. [DO Collaboration], Phys. Lett. B 517, 282 (2001) |arXiv:hep-ex/ 0106059 . 

[7] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 69, 052003 (2004). 

[8] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 65, 091102 (2002). 

[9] J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 (1999) |arXiv:hep-ph79 905386 . 
[10] J. M. Campbell and R. K. Ellis, Phys. Rev. D 62, 114012 (2000) ar Xiv:hep-ph/0006304 . 
[11] J. Campbell and R. K. Ellis, Phys. Rev. D 65, 113007 (2002) |arXiv:hep-ph/0202176| . ~ 
[12] S. Cortese and R. Petronzio, Phys. Lett. B 253, 494 (1991). 
[13] T. Stelzer and S. Willenbrock, Phys. Lett. B 357, 125 (1995). 



34 



[14] A. P. Heinson, A. S. Belyaev and E. E. Boos, Phys. Rev. D 56, 3114 (1997). 

[15] S. S. Willenbrock and D. A. Dicus, Phys. Rev. D 34, 155 (1986). 

[16] C. P. Yuan, Phys. Rev. D 41, 42 (1990). 

[17] R. K. Ellis and S. Parke, Phys. Rev. D 46, 3785 (1992). 

[18] D. O. Carlson and C. P. Yuan, Phys. Lett. B 306, 386 (1993). 

[19] T. M. P. Tait, Phys. Rev. D 61, 034001 (2000) |arXiv:hep-ph/9909352|. 

[20] A. Belyaev and E. Boos, Phys. Rev. D 63, 034012 (2001) |arXiv:he~p h/0003260 . 

[21] M. C. Smith and S. Willenbrock, Phys. Rev. D 54, 6696 (1996) |arXiv:hep-ph/9604223| . 

[22] G. Bordes and B. van Eijk, Nucl. Phys. B 435, 23 (1995). 

[23] T. Stelzer, Z. Sullivan and S. Willenbrock, Phys. Rev. D 56, 5919 (1997) 
|arXiv:hep-ph/9705398] . 

[24] B. W. Harris, E. Laenen, L. Phaf, Z. Sullivan and S. Weinzierl, Phys. Rev. D 66, 054024 

(2002) |arXiv:hep-ph/0207055| . 
[25] Z. Sullivan, |arXiv:hep-ph"704~0 8049 , 

[26] V. S. Fadin, V. A. Khoze and A. D. Martin, Phys. Lett. B 320, 141 (1994) 

|arXiv:hep-p h /9309234] . 
[27] V. S. Fadin, V. A. Khoze and A. D. Martin, Phys. Rev. D 49, 2247 (1994). 
[28] K. Melnikov and O. I. Yakovlev, Phys. Lett. B 324, 217 (1994) |arXiv:hep-ph/9302311| . 
[29] R. Pittau, Phys. Lett. B 386, 397 (1996) ||arXiv:hep-ph/9603265|. 
[30] C. Macesanu, Phys. Rev. D 65, 074036 (2002) |arXiv:hep-ph/0112142| . 
[31] R. K. Ellis, D. A. Ross and A. E. Terrano, Nucl. Phys. B 178 (1981) 421. 
[32] S. Catani and M. H. Seymour, Phys. Lett. B 378 (1996) 287 |hep-ph/9602277l . 
[33] S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratum-ibid. B 510 (1997) 

291] lhep-ph/9605323"] . 

[34] S. Catani, S. Dittmaier, M. H. Seymour and Z. Trocsanyi, Nucl. Phys. B 627, 189 (2002) 
|arXiv:hep-ph/0201036] . 

[35] Z. Nagy and Z. Trocsanyi, Phys. Rev. D 59, 014020 (1999) [Erratum-ibid. D 62, 099902 

(2000)] |arXiv:hep-ph/9806317| . 
[36] Z. Nagy, Phys. Rev. D 68, 094002 (2003) |arXiv:hep-ph/0307268| . 
[37] M. Jezabek and J. H. Kuhn, Nucl. Phys. B 314, 1 (1989). 
[38] T. Gottschalk, Phys. Rev. D 23, 56 (1981). 

[39] C. R. Schmidt, Phys. Rev. D 54, 3250 (1996) |arXiv:hep-ph/9504434| . 
[40] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 44, 189 (1972). 

[41] Z. Bern, A. De Freitas, L. J. Dixon and H. L. Wong, Phys. Rev. D 66, 085002 (2002) 

|arXiv:hep-ph/020227l1 . 
[42] G. Altarelli, R. K. Ellis and G. Martinelli, Nucl. Phys. B 157, 461 (1979). 
[43] Z. Nagy, private communication. 

[44] S. Dittmaier, Nucl. Phys. B 565 (2000) 69 |hep-ph/9904440] . 

[45] V. A. Khoze, W. J. Stirling and L. H. Orr, Nucl. Phys. B 378, 413 (1992). 

[46] P. Azzi el al. [CDF Collaboration] , |arXiv: hep-ex/ 0404010} 

[47] G. C. Blazey et al, arX iv:hep-ex/0005012| 

[48] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 

012 (2002) arXiv:hep-ph/0201195| . 
[49] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28, 455 (2003) 

|arXiv:hep-ph/0211080] . 



35 



